Anomalies of ac driven solitary waves with internal modes: 
Non-parametric resonances induced by parametric forces 



O 

o 

On 



O 



■4— > 

CO 

+-* 

a 

c 
o 
o 



> 

On 

m 

o 



a 



-a 
c 
o 
o 



X 
S3 



Niurka R. Quintero^ 

Grupo de Fisica No Lineal, Departamento de Fisica Aplicada I, Universidad de Sevilla, 
Facultad de Informdtica, Avenida Reina Mercedes s/n, 41012 Sevilla, Spain 

Angel Sanchez^ 

Grupo Interdisciplinar de Sistemas Complicados (GISC), Departamento de Matemdticas, 
Universidad Carlos III de Madrid, Avda. Universidad 30, 28911 Leganes, Madrid, Spain 

Franz G. Mertens § 
Physikalisches Institut, Universitdt Bayreuth, D-95440 Bayreuth, Germany 

(Dated: February 1, 2008) 

We study the dynamics of kinks in the 4 model subjected to a parametric ac force, both with 
and without damping, as a paradigm of solitary waves with internal modes. By using a collective 
coordinate approach, we find that the parametric force has a non-parametric effect on the kink 
motion. Specifically, we find that the internal mode leads to a resonance for frequencies of the 
parametric driving close to its own frequency, in which case the energy of the system grows as 
well as the width of the kink. These predictions of the collective coordinate theory are verified 
by numerical simulations of the full partial differential equation. We finally compare this kind 
of resonance with that obtained for non-parametric ac forces and conclude that the effect of ac 
drivings on solitary waves with internal modes is exactly the opposite of their character in the 
partial differential equation. 



I. INTRODUCTION 

During the last three decades of the 20th century, 
we have witnessed a great deal of effort devoted to the 
study of solitons and related nonlinear coherent excita- 
tions |y, Q, ||. All that work notwithstanding, the analy- 
sis of systems generally regarded as paradigms of soliton 
phenomenology still yields unknown, surprising features. 
One such case is the dynamics of topological solitons or 
kinks in nonlinear Klein-Gordon systems subjected to ac 
forces. After several papers that led to contradictory con- 
clusions during the last decade, the problem for the sine- 
Gordon equation was finally solved in 1998 M. This was 
the simplest possible scenario for studying the effects of 
ac forces on topological solitons, in so far as sine-Gordon 
kinks do not have any internal mode and behave as rigid 
objects M. However, this is a nongeneric situation, be- 
cause many solitary waves do possess internal degrees of 
freedom || g, |7J , and even systems like the sine-Gordon 
model can acquire internal modes due to the influence of 
perturbations, such as discreteness [pL 0, M. Therefore, 
it was necessary to study the behavior of solitons with in- 
ternal modes, and subsequently it was found that pure, 
non-parametric ac forces induced resonances of paramet- 
ric character: If the frequency of the internal mode is 
denoted by fli, the resonance appears when the driving 
frequency is fij/2 111). This anomalous behavior was 
shown, taking the (^equation as a typical example, to 
arise from the fact that the driving couples indirectly to 
the internal mode, which in turn leads to a parametric 



influence on its evolution. 

In view of the above results, it is only natural to pose 
the following question: If non-parametric drivings act 
parametrically on <j> kinks, what is the effect of a para- 
metric driving? To our knowledge, this issue was first 
addressed in [Ol , although the main point of that work 
was to show how kinks can be effectively anihilated by a 
fast parametric driving, and no attention was paid to the 
kink dynamics. More recently, the problem was consid- 
ered with both parametric and non-parametric forces act- 
ing simultaneously in the system (13], but the approach 
employed by those researchers was not appropriate and 
led to incorrect results W%- Therefore, the dynamics 
of kinks in the </> 4 model subjected to parametric driv- 
ings is largely unexplored, and this is another reason 
why we concern ourselves with this problem. This pa- 
per presents our conclusions on this question within the 
following scheme: Section II is devoted to our analyti- 
cal approach; Section III contains the simulations of the 
full partial differential equation that verify our analytical 
predictions, and Section IV summarizes our main results. 



II. RESONANCES INDUCED BY 
PARAMETRIC PERIODIC FORCES 



As we are interested in the behavior of solitary waves 
with internal modes, we choose as our working example 
the c/> 4 equation, which is well known to be representative 
of the generic behavior of those solitary waves. Therefore, 



we begin by considering the perturbed problem given by 

<ht-4>x* = —jT-P>l>t + f(t,<i>), (i) 

where U{4>) = (1 - 2 ) 2 / 4 is the 4> A potential, f(t,<f>) = 
esm(5t + §o)4> represents the parametric ac force, e, (5, 
and Sq are the amplitude, frequency, and phase of the 
external periodic force, respectively, and (3 is the damp- 
ing coefficient. 

In order to obtain analytical results, we resort to the 
well-known collective coordinate method 0, Il5), which 
will give us definite predictions on the behavior of the 
system. However, these being approximate results, we 
will have to check them by direct numerical simulations 
of Eq. (Q); those will be the subject of the following sec- 
tion. To apply the collective coordinate approach to our 
problem, we choose to employ the so-called generalized 
traveling wave Ansatz (T6( or, equivalently, to use the 
variation of the momentum and the energy (the equiva- 
lence of both techniques is proven in the second paper in 
PI). These two quantities are given by 



P(t) 
E(t) 



dxcj> x (f>t, 
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In order to evaluate those quantities, we consider the Rice 
Ansatz llTll and assume that 



We now have to evaluate the above quantities. In this 
respect, it is interesting to compare the equations we have 
obtained with those found for the case of an ac non- 
parametric force, i.e. f(t,(f>) = f{t). In that situation 
we obtain that F ex = —qf(t) and K = 0, and hence the 
external force acts directly on the translational mode, 
whereas the internal mode is excited indirectly due to 
the coupling between these two modes. Conversely, as 
we have seen, in the parametrically driven problem, F ex 
vanishes, while K = — esin(5t + <5 ). Therefore, Eqs. (g)- 
(|6|) become 



f = ~ m - 
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We thus see that, when introduced in a parametric man- 
ner, the ac force acts directly on the internal mode and 
not via the coupling between the translational and inter- 
nal modes. This result shows the very different role of 
both drivings at the collective coordinate level. 

In order to deal with these equations, we choose 
X(0) = u(0), 1(0) = l s = Wl-« 2 (0) and 1(0) as the 
initial conditions for Eqs. ( |l0|) and (|ll|). The equation 
for the momentum is trivial and can be solved exactly, 
yielding 



(x, t) = 0o 



X(t) 



l(t) 



(4) 



where <po(x) is the static kink solution of the unperturbed 
(j) A equation, i.e., of Eq. (nh with (3 = e = zero, centered 
at X(t) and with width l(t). With this choice, and fol- 
lowing the standard procedure [|l Q5[, we find that the 
evolution of the momentum, P(t) = MqIqX /l(t), and the 
width of the kink are given by 
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where q = 2, Mo = 4/(3fo)) and Iq = \/2 represent the 
topological charge, the mass, and the width of the un- 
perturbed kink, respectively, and a = (ir 2 — 6)/12 is a 
constant. 



P(t) = P(0)e 



-0t 



(12) 



We note that for (3 ^ 0, the momentum of the kink goes 
to zero and after some transient time, t » 1/(3, it ef- 
fectively vanishes. In the same way, if we start from zero 
initial velocity, P(t) — 0, i.e., the center of the kink will 
not move. 

Let us now turn to Eq. (pT|). As a first step, we in- 
troduce a change of variables, proposed in |[I8|, given 

111) into 
which 



by l(t) = g 2 (t). This change transforms Eq L 



an Ermakov-type equation (or Pinney-type) 
reads 
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where fl = I/^/oHq = 1.2452 is equal to the Rice fre- 
quency |TiJ Qr = 1/ ' \fal s in the case when the kink 
initially is at rest. It can be shown (see |17 and [p0| ) 



that when there is not any perturbation in the system 
(H)-(O), l(t) oscillates with a frequency flu if we start 
from any ^(0) ^ l s or 1(0) ^ 0. Therefore, flu is the 
characteristic frequency of Eq. (O) and, since it agrees 



within 1.7% with fli = ^/ZJ2 = 1.2247, we expect that if 
we find a resonance related with fin in Eq. (13J) [or Eq. 
(O)] we should find a similar phenomenon in the full sys- 
tem, Eq. (|l|), associated to the frequency of the internal 
mode Qi. 

We now proceed to solve the Eq. ( |l3|) analytically for 
the undamped case, i.e., when (3 = 0. As we have men- 
tioned, in this case P(t) — P(0) and Eq. ( |l3|) becomes 
the following Pinney-type equation p9[ 



2 
whose solution is 



sin (ft + Sq) 
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(16) 



where v\ (t) and v 2 (t) are two independent solutions of the 
linear part of Eq. (Rq) and W — V\V 2 — v 2 V\ is the Wron- 
skian. W(t) is actually a constant, and can be calculated 
from the initial conditions for Vi (i — 1,2), i>i(0) = \/T s , 
Vi(0) = Z(0)/(2v / ?I), v 2 (0) = and 1)2(0) a nonzero con- 
stant. 

If one denotes t = (ft + <S + 7r/2)/2, after some alge- 
braic manipulations we arrive at the following Mathicu 
equation for the Vi functions: 



5 



[a - 26 cos(2r)]vi = 0, 
aM W 2 Mo7o 



X) ■ il ' 1 



where prime denotes the derivative with respect to r, 
and 70 = 1/yl — t>(0) 2 . Notice that the initial condi- 
tions, when t = To = (<5o + 7r/2)/2 for Uj(r), become 
«i(r ) = V£, «i(r ) = i(0)/(Sy/Q, v 2 (t q ) = and 
w^(r ) = 2v 2 (0)/6. The solution of Eq. @ (see j|l)) 



for «i(r) and i>2(t) can be expressed as a linear super- 
position of the two Mathieu functions ce^ and se^ with a 
non-integer index v, i.e., 



Vi(r) = Aice v (t, < 
where 



Bi se u (r, ( 



1,2, 



Ai = 



A.. 






(18) 



(19) 



and 



A = ce„ (t ,9) se^ (t 0i 9) - ce' v (r ,6) se„ (t q ,9), 

A^i = Vi(ro) se^ (r , 9) - w-(t ) se„ (r , 0), 

As ; = v[(tq) ce„ (r , 0) - Ui(r ) ce'^ (r , 0), 

with the constraint (characteristic curve for Mathicu 
functions) 



2(i/ 2 - 1) 



0(9 4 ). 



(20) 



From Eqs. (^6|), ( Jlq ) and ([19]), and taking into account 
that r = (St + 5o +7r/2)/2 we obtain that the kink width 
l(t) is given by 



l(t)=g 2 =v 2 (t) + 



1 



4aIF 2 



v 2 2 (t), 



(21) 



where 



i7i(i) = A, ce„ (ft/2 + (5 /2 + tt/4, 
+ B s ; se„ (St/2 + S /2 + tt/4, 

w = -Vw> 2 (o), 



i = l,2;(22) 

(23) 



and the characteristic curve, Eq. (|20|), for our initial pa- 
rameters can be written up to order e 2 as 



V HI 



1/ I 4(^-l)M 2 7o 41 



+ 0(e 4 ). 



(24) 



Notice that when j/ = m + p/s is rational, with m an 
integer number and p/s a rational fraction (0 < p/s < 1), 
vi(t) and i>2(£) are 2-7TS periodic functions, if p is odd, and 
7rs-periodic functions, if p is even, whereas for irrational 
v both functions will be non-periodic although bounded 
pl[ . For instance, if we take in Eq. (11) e = 0.01, 
do = vr/2, u(0) = 0, 1(0) = l s and S = 0.94, we can 
see in Fig. [j] that l(t) is a bounded function. Moreover, 
taking the Mathieu functions that appear in Eq. ( |22| ) 
up to terms of order 9 and substituting the approximate 
expression of Eq. (122) in Eq. (Efi) one can show that 
the resulting expression for l(t) involves basically the fre- 
quencies FIr and S (see Fig. nb. We have also plotted 
in Fig. H the evolution of the energy, Eq. (||), taking 8 
not close to £Ir. For this choice of parameters, E(t) is 
also a bounded function and its spectrum for the above 
parameters involves chiefly three frequencies: \8 — £Ir\, 
28 and 8 + flu. 

As we have shown, Eq. (pil) represents the charac- 
teristic curve corresponding to stable (either periodic or 
bounded) solutions of the Mathieu equation (17). How- 
ever, if we try to find a periodic solution when 8 w flu, 
we obtain that the integer Mathicu functions sei and cei 
are two independent solutions of two different Mathicu 
equations since they are related to different characteristic 
curves: 



a = 1 — 1 

0=1 + ' 



04 

03 



1536 

fl4 



8 64 1536 



(25) 
(26) 



respectively. These characteristic curves separate the un- 
stable and stable regions of Eq. (|17|). In the unstable 
regions the solution of this equation, as well as l(t) and 
E(t), grow with t (see Fig. g). This is hence the hint 
on the existence of a resonance phenomenon, and, as we 
already pointed out, we expect that it will manifest itself 
in the full partial differential equation with an increment 
of the width and the energy of the kink, when 8 is close 



to Qi. This is the aim of the next section, in which we we 
will check by computing the solution of (|l|) whether there 
is some resonance at 8 w fij. In addition, we will analyze 
numerically the effect of damping in our resonance pic- 
ture, as we have not been able to obtain any analytical 
approximate result in that case. 



III. SIMULATION RESULTS 



In this section we carry out numerical simulations of 
Eq. (|lj) by using the Strauss- Vazquez scheme |22] , taking 
the length of the system L = 100, At = 0.01, Ax = 0.1, 
the final time t/ = 200 and free boundary conditions. We 
fix also the amplitude e = 0.01 and the phase So = n/2 
of the ac force. 

In order to show the evolution of l(t) and E(t) in the 
PDE (|l|) we distinguish two cases: In the first case we 
choose 8 far away from f2j, whereas in the second one 
8 w f2j. In Fig. we plot the width and the energy of the 
kink as functions of t for 8 — 0.94, i.e., the parameters of 
the ac force are the same as in Fig. pi If we compare Figs. 
HI and [| we see that l(t) and E(t) have the same behav- 
ior both on the collective coordinate level and in the full 
system, Eq. (Q). Naturally, while in the spectrum of /(t) 
obtained by the collective coordinate approach there ap- 
pear the frequencies 8 and £Ir, in the DFT of l(t) for the 
PDE we obtain 6 and fij as expected. However, in the 
Fourier transform of E(t) for the full system we find not 
only frequencies which involve 8 and f2i (see Fig. ||), but 
also an additional frequency, lo^ = 0.4715 ~ ujq — 8, re- 
lated with the frequency of the lowest phonon (u = \/2) 
of the (j) A system. This means that the internal mode 
and the phonons can both appear when the <j> system 
is driven parametrically with an ac force, but since fli 
and too are separated we expect to be able to excite the 
internal mode to a larger extent than phonon modes if 
the frequency of the ac force is closer to fii than to ujq. 

Moving now to the second case, i.e., drivings with fre- 
quency 8 close to f2j, Fig. E shows that in this situation 
the width of the kink oscillates with a large amplitude 
while the energy grows, although not monotonously. It 
is important to stress at this point that, whereas the 
collective coordinate analysis predicts that l(i) and E{t) 
should increase indefinitely at resonance, this is not the 
case in the original system. The reason is that, in our 
perturbative method, we have not taken into account the 
phonons, which are present in the system, and therefore 
part of the energy input from the ac force goes to the 
phonons. In our simulations we have also seen that the 
kink does not move, even if we are at resonance, in agree- 
ment with our collective coordinate theory. 

In order to present in a more evident manner the res- 
onance, we have plotted in Fig. H the values of the mean 
normalized width and the mean energy as functions of 
8. We compute these magnitudes as time averages of 
l(t)/lo and E(t) for times 100 < t < 200, as to avoid pos- 



sible transient effects (which are mostly relevant in the 
damped case discussed below). We see that on the col- 
lective coordinate level (l(t)/lo) and (E(t)) have a maxi- 
mum at 8 = Qr] correspondingly, in the original system, 
Eq. ([[]), the resonance takes place at 6 = 1.18, close to 
fii = 1.2247. We have also found that for 1.3 < 8 < uj 
(not shown), the energy is an increasing function of 8 
whereas the width of the kink is practically constant. 
This implies that the increment of the energy is due to 
excitation of the phonons, which becomes more efficient 
as 8 — » wo, and is not related with the internal mode. 

Coming now to the damped case, as we announced in 
the preceding section, we have not obtained any analyt- 
ical result. For this reason, in this section we carry out 
numerical simulations of Eq. (gj) and also integrate nu- 
merically the collective coordinate prediction, Eq. (pi]), 
when /3 ^ 0. In Fig. J6| we plot the evolution of the energy 
at resonance for (3 — 0.05, as given by the collective coor- 
dinates and as obtained from the numerical simulations 
of Eq. (Q). For each case, we have chosen a frequency 
near to fi/j (for the collective coordinates) and fii [for 
Eq. p)]. To complete the study of the damped case, we 
present the result for the resonance as we did for the un- 
damped case, by plotting the mean width and energy vs 
8 for j3 — 0.05 (see Fig. [7]). We can observe that, at reso- 
nance, these functions also have maximums, although in 
the presence of damping they are less peaked than in the 
undamped case: (l(t)/lo) increases around a 7% and the 
energy around a 10%. Interestingly, the resonance for the 
full system is now much closer to the internal mode than 
in the undamped case. We believe that the fact that the 
damping suppresses quite effectively the phonons is the 
reason why the damped system follows more closely the 
collective coordinate prediction, which indeed neglects all 
phonon contributions. 



IV. CONCLUSIONS 

In the present paper we have studed the effect of a 
parametric ac force on the dynamics of </> 4 kinks with 
and without damping. By using a collective coordinate 
approach, we have found resonances related with the ex- 
citation of the internal mode when the frequency of the 
ac force is close to the Rice frequency flu sw £V We have 
verified numerically this prediction by computing the so- 
lution of the Eq. (pi), its corresponding energy and the 
width of the kink. When this kind of resonance occurs, 
the energy and the width of the system increase, while 
for other values of the frequency (except for smaller, sec- 
ondary resonances) they are bounded (in some cases pe- 
riodic) functions. Concerning the mobility of the kink, it 
is very important to point out, at least for the time con- 
sidered here, that all the input energy goes to the internal 
and phonon modes and for this reason the kink does not 
move (neither near nor away from the resonance) . 

Although this work has focused on the 4 kink, ex- 



perience shows that it is very likely that similar phe- 
nomenologies will be observed in the dynamics of other 
solitary waves with internal modes. In this respect, it 
is important to compare the above results to those ob- 
tained for non-parametric ac drivings [rifl in order to give 
a thorough, coherent picture of the problem. The col- 
lective coordinate analysis, confirmed by the direct nu- 
merical simulations of the full system, shows that when 
we drive the <j) A system parametrically the force acts di- 
rectly and non-parametrically on the width of the kink, 
i.e., the internal mode is excited and the phenomenon of 
resonance takes place and can be observed by monitoring 
l(t) and E(t). On the contrary, non-parametric drivings 
affect the internal mode only indirectly through the cou- 
pling to the translation mode, which in turn leads to a 
parametric coupling. We thus see that the internal mode 
dynamics is responsible for a highly non-trivial behavior: 
non-parametric drivings on the full system induce para- 
metrical resonances in the kink dynamics, and vice versa. 
As a final remark, we want to point out that our results 
indicate a way of exciting internal modes without affect- 
ing the translation mode, which can mask the effects one 
is interested in. As we have seen, the parametric ac force 
affects only the internal mode, and the kink remains at 
its original position (or moving uniformly with the same 
speed if nonzero) at all times. This is not possible with a 
non-parametric driving and might be a mechanism of in- 
terest in situations where the role of internal modes must 
be elucidated. 
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FIG. 1: Collective coordinates: Evolution of l(i) and E(i) 
when 5 = 0.94 for the undamped case, j3 = 0. Points: 
Analytical solution [see Eqs. (|lq )-(|2^)1. Solid lines: Nu- 
merical integration of Eq. (|ll|) for j3 = 0, e = 0.01, 
6 = vr/2, u(0) = 0, X(0) = 0. The DFT (not shown) 
of l(t) yields two frequencies: u)\ = 0.9434 rs 5 and 
W2 = 1.25789 fs f2#, whereas for £'(i) it gives three 
frequencies: Wi = 0.3144 fs \S - Sl R \ = 0.3052, uj 2 = 
1.8868 « 25 = 1.88 and w 3 = 2.20131 ^§ + Q, R = 2.1852 
. The maximum of the DFT appears at io\ . 
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FIG. 2: Collective coordinates: Evolution of l(i) and E(t) 
when 5 = 1.24 ps fi# for the undamped case, (3 = 0. 
Points: Analytical solution [see Eqs^ ©- @]. Solid 
lines: Numerical integration of Eq. 
as in the previous figure. 
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FIG. 3: Evolution of l(t) and E(t) as obtained from 
numerical simulations of Eq. (|l|), when 5 = 0.94 with 
e = 0.01, S = tt/2, w(0) = 0, X(0) = 0. The DFT (not 
shown) of l(t) yields two frequencies: u>\ = 0.9429 sw 
8 and lo^ = 1.2258 w Oi, whereas for E'(i) gives four 
frequencies: w x = 0.2887 w |<5 - fij| = 0.2847, w 2 = 
1.8859 w 2(5 = 1.88, cj 3 = 2.3574 w S + n t = 2.1647 and 
uj 4 = 0.4715 «w -i = 0.4749. 
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FIG. 4: Evolution of l(t) and E(t) as obtained from 
numerical simulations of Eq. (fiT), when 5 = 1.176. Other 
parameters are the same as in the previous figure. 
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FIG. 5: (l(t)/lo) (points connected by solid line) and 
(E(t)) (dashed line) vs 5. Upper panel: Collective co- 
ordinate approach. The resonance is at flu = 1.2452. 
Lower panel: Results from the full system, Eq. (|l|). The 
resonance is at S = 1.18, close to fL = 1.2247. 
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FIG. 6: Damped case, /3 = 0.05: Evolution of the energy 
according to the collective coordinate approach, obtained 
by numerical integration of Eq. (|ll|), and from the full 
system, by numerical simulation of Eq. p), when 5 = 
1.245 w £Ir and S = 1.22 « O^ respectively. 
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FIG. 7: Damped case, (3 = 0.05: (l(t)/l ) (points con- 
nected by solid line) and (E(t)) (dashed line) vs 5. Up- 
per panel: Collective coordinate approach. The width 
and the energy take their maximum values at S = 1.246 
and 6 = 1.245 respectively. Lower panel: Results from 
the full system, Eq. ([[]) . The effective width has its max- 
imum at 5 = 1.22 and the maximum of the energy occurs 
at S = 1.21, both of them close to fl t = 1.2247. 



